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Abstract 

We present a deterministic algorithm that given any directed graph 
on n vertices computes the parity of its number of Hamiltonian cycles 
in 0(1. 618") time and polynomial space. For bipartite graphs, we give 
a 1.5" poly(w) expected time algorithm. 

Our algorithms are based on a new combinatorial formula for the num- 
ber of Hamiltonian cycles modulo a positive integer. 



1 Introduction 

It is known since the 1960s that Hamiltonian cycles in an n-vertex graph can 
be detected and counted in 0{2"n^) time |1 8|. In an influential survey, Woeg- 
inger |11] asked if this could be significantly improved. Recently, two results 
using different techniques have met this challenge: 

1. a randomized 0(1.657") time algorithm for the decision problem in undi- 
rected graphs (2), 

2. a randomized 0(1.888") time algorithm for the decision problem in di- 
rected bipartite graphs f5\. 

For the general directed graph case, no such algorithm has appeared yet. 

Inspired by the fact that the algorithm for the bipartite case in |5| actually 
computes the parity of the number of Hamiltonian cycles, we show how to 
compute this value for any directed graph in 0(1.618") time. 

While this does not seem to shed any direct light on the decision problem 
in directed graphs, it does however ask a related interesting question: Could 
it really be easier to solve a ®P-complete problem than its decision counter- 
part? Even the known algorithm for the undirected graph decision case [2] is 
slower than our present algorithm. We note in passing that for another well- 
known problem. Set Cover, current evidence points in the opposite direction: 
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It is known that a fast algorithm for computing the parity of the number of set 
covers would disprove the Strong Exponential Time Hypothesis [4J, whereas 
no such connections are known for the decision variant. 

1.1 Results 

Let denote the set of (directed) Hamiltonian cycles of a directed input graph 
G = (V, E) on n = \V\ vertices. 

Our main result is that the parity ®^ of the number of directed Hamilto- 
nian cycles can be computed much faster than in 2" time. 

Theorem 1. We can compute (B,^ for a general n-vertex digraph in time within a 
polynomial factor of the nth Fibonacci number F„ e 0(1.618") and polynomial space. 

For the restricted family of bipartite graphs, we can show a stronger boimd. 
This result is also somewhat easier to prove. 

Theorem 2. We can compute for a bipartite n-vertex digraph in expected time 
within a polynomial factor ofO{1.5") and polynomial space. 

Both results rely on a new characterization of the number of Hamiltonian 
cycles in terms of a local property. 

Theorem 3. For a vertex subset X, let di,(X) denote the number of directed edges 
from vertex vto a vertex in X. For integer K >2, 

m = ^ E i--^)^''''^^(Udz{V\Z))Yl(Yld,{Y,)) (modK), 
z,Yi,...,Yk \zez / k=i Vgy^. / 

where the sum is over all {K+ 1) -partitions ofV. 

Our algorithms use this result for K = 2, but we state (and prove) it for 
general K. 

1.2 Related results 

There are nontrivial examples of hard counting problems where computing the 
parity is easy. In fact an important example is given by a closely related prob- 
lem: computing the number ft*^ of disjoint cycle covers in a directed graph. 
This is a hard counting problem, equivalent to computing the permanent of 
the adjacency matrix. On the other hand, the parity ffi*^ equals the parity of 
the determinant and is therefore computable in polynomial time by Gaussian 
elimination. However, the problem studied in the present paper, re- 
mains complete for the complexity class ©P, even for very restricted classes 
of graphs [TO]. 

Counting the number of Hamiltonian cycles in an «-vertex directed graph 
can be done in o(2") time (Sj. However, the improvement over the classic al- 
gorithms (TJISI is relatively small (within a factor of exp(0(-\/«/ log«))). The 
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parity problem appears to be easier. First, for undirected graphs it is known 
that when all vertex degrees are odd, there is an even number of Hamilto- 
nian cycles passing through each (directed) edge |9|. Second, a recent paper of 
Cygan, Kratsch, and Nederlof |5| includes a deterministic 0(1.888") time, ex- 
ponential space algorithm computing the parity of the number of Hamiltonian 
cycles for undirected and directed bipartite graphs. Our present algorithm is 
faster, uses only polynomial space, and works for unrestricted directed graphs 
as well. On the other hand, the algorithm from |5| works for weighted graphs. 
This is a crucial property because it allows the use of the Isolation lemma to 
construct a (randomized) decision algorithm. Our constructions do not seem 
to allow this extension. 

The parity and decision problems can be compared on the promise prob- 
lem of detecting a unique Hamiltonian cycle. The result of [2] provides a 
randomized algorithm for undirected graphs; our result also applies to di- 
rected graphs, and is slightly faster on undirected graphs (0(1.618") versus 
0(1.657")). 

1.3 Open questions 

Our paper accentuates the fact that the exponential time complexity of count- 
ing or deciding the Hamiltonian cycles in a directed graph is not very well im- 
derstood. While there are several examples of ®P-complete problems whose 
decision analogue is computationally easier (e.g., ®2-Satisfiability), examples 
of the converse are not known to the authors. 

Under the Exponential Time Hypothesis, the decision problem does not 
allow exp(en)-time algorithms for arbitrary e > 0, but there are no arguments 
for or against an 0(1.999") -time algorithm for the decision problem. 

For comparison, it is now known that the parity of the number of set covers 
of a given set system on n elements requires 0(2(^^'^'") time under the Strong 
Exponential Time Hypothesis 14J. 

1.4 Overview 

We first describe the top level algorithm in section |2l based on theorem |3l In 
sections |3] and |4]we address the specific parts of the algorithm for the bipartite 
and the general input case, respectively. Finally, in section|5]we give the proof 
of theorem m 

2 Algorithm 

By specializing the expression from theorem|3]to K = 2, we get 

®-^=2 E (Ud-i^\^))(Udyiyi))(Udym) (mod 2) . 
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In particular, we removed the factor ( — because —1 = +1 (mod 2). 

We introduce X = Yi, Y = Y2 and let Z denote the vertex complement 
Z = y \ Z. Finally, we separate the expression into the contributions of X and 
Y U Z, respectively: 

(BJi^ = Y^f{X)Y[d,{X) (mod 2), (1) 

where 

fix) = ij:(uMy))udz(z). (2) 

Y,z \ijeY / zez 

Still, the sum is over all tripartitions X, Y, Z, of the vertices and thus has 3" 
terms, so it does not in itself serve as a fast algorithm. 

There are two parts to our algorithm. First, for many X, already the factor 
Uxex dxiX) in ^ will be zero modulo 2. Hence, these subsets contribute noth- 
ing and we will avoid explicitly listing them. Second, for any X for which it is 
not zero, the value of /(X) modulo 2 can be computed in polynomial time in 
|y| via linear algebra. 

On the top level our algorithm is simply an evaluation of and 

Algorithm P (Parity.) Given a directed graph G = {V, £), compute © J^. 
PI [Initialise.] Set s = 0. 

P2 [Locate.] List every X C V such that U^^x^xiX) is odd. 

P3 [Contribute.] Compute /(X) (mod 2) for every such X and add it to s. 

P4 [Report.] Return s (mod 2). 

Note that steps F2 and P3 can be interleaved in order to avoid storing all X. 
We proceed by explaining how to execute steps P2 and P3 efficiently. 

2.1 Step P2: Locating Subsets that Contribute 

We say that a subset X C V is contributing if ]l^^gx dxiX) is odd. Our aim is to 
generate the contributing subsets. We identify V with {1, . . . , n} and introduce 
the indicator variables Xi,...,x„ e {0, 1} for X C V with a:, = 1 if and only if 
i E X. Then the constraint that d, (X) is odd can be expressed in terms of the 
adjacency matrix A of G as 

n 

UijXj = 1 (mod 2) for all z e X, 

equivalently, 

Xi dijXj^ = Xj (mod 2) for all z = 1, . . . , n . 
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We will view these constraints as a system of n quadratic equations over 
GF(2) in the variables xi, . . . , 

+ ai2X2 + ■ ■ ■ + ai„Xn) = Xi 
^2(^21^1 + ^22^2 + ■ ■ ■ + Clln^n) = ^2 

X„{a„iXi + fl„2^2+ ■ ■ ■ + UnnXn) = ■ 

This system can be succinctly expressed as a matrix equation over GF(2): 

X o Ax = X (3) 

where x = {xi, . . . ,Xn) and the operator o denotes the coordinate-wise (or 
Hadamard, or Schur) product. 

2.2 Random self -loops 

A priori, listing the solutions to ^ requires 2" steps. In fact, already the size of 
the solution set is easily seen to be of order 2" for some graphs. For instance, 
every vector x with an odd number of Is solves (|3) if A is the all-Is matrix 
corresponding to a directed clique with self-loops. 

This motivates the main algorithmic insight in our construction. Because 
no Hamiltonian cycle can use a self-loop, we can add self-loops to the vertices 
G (or remove them) without changing jyf. Algebraically, manipulation of self- 
loops corresponds to flipping the diagonal entries of A. 

We show that adding these loops at random reduces the expected solution 
size for all graphs. 

Lemma 1. Let A be an n x n matrix with 0, 1-entries. Choose the diagonal entries 
ail,. • ■ ,cinn G {0/1} uniformly and independently at random. Then the expected 
number of solutions to ^ is 1.5". 

Proof. Consider the fth entry of x. If x/ = then the z'th equation in Q is = 
and trivally satisfied. If X/ = 1 then the fth equation in (|3]l is satisfied if 

n n 

1 = E = «ii ■ 1 + E '^iFi ' 
j=i j=i 

which happens with probability j because fl;; is uniformly distributed. By in- 
dependence of the choices of an, a vector x with k Is satisfies all the equations 
with probability 2"*^. By linearity of expectation, the expected number of solu- 
tions is 

|:(:)^-'Mi+ir, 

using the binomial theorem. □ 
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In section|3]we show that these 1.5" solutions can be efficiently listed if G is 
bipartite. In section|4]we show that the solutions can be listed in time 0(1.618") 
for general graphs. 



2.3 Gaussian elimination 

We assume that the following algorithm for solving nonhomogeneous systems 
of linear equations is well known. We repeat it here only to recall that the 
solution set can be described in terms of the basis of translated vector space 
(the system's null space). In particular, if a linear equation system has a solution 
at all, the number of solutions equals 2'^, where d is the dimension of the null 
space. 

Algorithm G {Gaussian elimination.) Given a nonhomogeneous system Ax = c 
of linear equations. If the system has no solution, output "no". Otherwise, 
output a vector v and d linearly independent vectors fci, . . . , b^j such that the set 
of vectors x with Ax = c equals {v + span(foi, . . . ,b^)}. 

Gaussian elimination runs in polynomial time and works over finite fields. 



2.4 Step P3: Computing the Contributions 

We return to the value /(X) in ||2) and show how to compute it in polynomial 
time, at least modulo 2. 

It will be convenient to remove the factor j . To this end, consider a partial 
order -< on the subsets of V such that for disjoint X, Y C V, not both empty, 
either X -< Y or Y ^ X. Then, from 0, 

Y,z\ijeY / zez 



Y,Z \jeY / Z6Z 

E (Udyiy)) n rfz(XUY), (4) 



Ycy\X ^i/GY / ziXUY 

X-<Y 

where Y, Z always form a disjoint partition of V \ X. To see the second step, 
note that we never get a contribution for Z = V because the term Ylzez dziZ) 
vanishes in this case. Hence, X U Y is non-empty for any contributing Z, and 
because of the symmetry in the formula with respect to X and Y, each Z that 
contributes to the sum will be counted an even number of times. 

The salient feature of @ is that modulo 2, we can now view /(X) as the 
number of solutions to an equation system over GF(2). 
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We define the partial order -< on the subsets of V as follows. Identify V 
with {!,...,«} and set X ^ Y if min X < min Y, with the usual convention 
minO = 00. (By the previous paragraph, the case X = Y = never arises.) 

We will set up a system of linear equations in 1/1, . . . , so that /(X) equals 
the number of solutions (mod 2). We use the obvious correspondence 



1, if/eY; 
0, if/^Y. 



First, since Y belongs to the complement of X we introduce the equations 

y, = f or / e X . (5) 

Second, to ensure X ^ Y we introduce the equations 

y, = 0, for z = 1, . . . , min X . (6) 

Finally, a subset Y with X -< Y contributes to /(X) if all product terms in 
(Sll are 1, so we need 

rf/(Y)=l (mod 2), iffeY 
di{X U Y) = 1 (mod 2), if f ^ Y 

for each i ^ X. Since d/(X U Y) = clj{X) + rf,(Y) we can rewrite the second 
constraint to 

d,(Y) = 1 + di{X) (mod 2), if ; ^ Y . 

In terms of the variables yi, . . .,y« and the adjacency matrix A of G we have 
the equations 

1 , if y,- = 1 ; 

l + di{X), ify,=0; 



for each i ^ X. (Note that d/(X) is a fixed value depending only on X.) The 
right hand side simplifies to 1 + d, (X) + d, (X)y„ so the resulting equations are 

n 

di{X)yi + E = 1 + di{X), for i^X. (7) 

In summary, the solutions to the equations in ||5]l, l|6ll, and ^ describe ex- 
actly the subsets Y C V \ X with X ~<Y that contribute 1 to /(X). The number 
of such solutions is odd exactly if the system of equations has a unique solu- 
tion. This is determined by Gaussian elimination in time pol5momial in n. 



3 Bipartite Input Graphs 

We first describe an algorithm for the step P2 for the case that the input graph 
is bipartite. We can safely assume that the input bipartite graph is balanced on 
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two equal sized vertex sets V-[ U 
nian cycles.) 
Let 

B = 



V. (Otherwise, there can be no Hamilto- 





B21 



B12 




denote the adjacency matrix of the input graph. Choose the random values 
ri,...,r„ E {0,1} uniformly and independently and construct the diagonal 
matrix R = diag(ri, . . . ,r„). Set A = B + R. The equation system Q becomes 



equivalently. 



xo (B + R)x = xo 



Xi O B12X2 
X2 O B21X1 



B12 
B21 



X + X o Rx = X , 



xi — xi o RiXi 

X2 — X2 O R2X2 / 



where x = (xi,X2) and R = diag(Ri,R2)- Iri this formulation it is apparent 
that for fixed xi, the equation system (in the variables X2) is linear; note in 
particular that X2 o R2X2 = K2X2 because R2 is diagonal and multiplication in 
GF(2) is idempotent. Thus, the system can be solved by Gaussian elunination 
in polynomial time. 

Algorithm B {Bipartite graphs.) Given a bipartite directed graph G = {V,E) 
with n = \V\, lists the vertex subsets X Q V for which di-(X) is odd for all 
X e X. 

Bl [Add random self loops.] Choose rj, . . . , r„ G {0, 1} uniformly at random 
and add a self -loop at vertex if r; = 1. 

B2 [Initialise.] Set xi = ■ ■ ■ = x„/2 = 0. 

B3 [Solve.] Solve the system x o Ax = x in the free variables x„/2+i/ ■ ■ ■ ,Xn 
using algorithm G. If it returns "none" proceed to B5. 

B4 [Report.] (Algorithm G returned the solution set as v,bi, . . . , b^.) Output 
the vectors v + ocibi + ■ ■ ■ + ocdbd for all 2** choices of ai, ■ ■ ■ , e {0, 1}. 

B5 [Next.] If all choices of (xi, . . . , x„/2) have been inspected, terminate. Oth- 
erwise generate the next choice and return to B3. 

The number of choices of xi is 2"^^, so step B3 takes total time 2"^^poly{n). 
This is dominated by the overall running time of step B4, which is linear in the 
number of of contributing subsets X. These are 1.5" in expectation according 
to lemma[TJ Hence the runtime is 1.5" poly(M) in expectation. This finishes the 
proof of theorem|2l 
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4 General Input Graphs 



For general graphs, we do not know how to list the solutions S in step P3 to 
(O in time proportional to the solution set. Instead, we will efficiently list a 
superset S' of the solution set of size 0(1.618"). We can then examine every 
X S S' to see that it solves (|3j. 

The superset S' of candidate solutions to l|3j is defined as follows. Let x' 
denote a prefix vector of fixed values x' = (xj, . . . , x„_j() E {0, Ij"^^ and intro- 
duce the variables . . . , x„. Set x = (xj, . . . , x„). Consider the equation 
system 

(x',0^)oAx= (x',Ofc). (8) 

For fixed x', this is nonhomogeneous set of linear equations in the variables 
x„_]^^l, . . . ,x„. Every prefix of a solution to l|3]l is a solution to 

Lemma 2. Let x = {xi,...,x„) e {0,1}" and x' = [xi, ... ,x„_{) for some k e 
{0, 1, ... , n}. Ifx satisfies x o Ax = x then (x', 0,^) o Ax = (x', Oj^). 

Proof. For 1 < i < n — k, the z'th equation of both x o Ax = x and (x', 0,^) o Ax = 
(x',Ofc) is 

n 

■^i 5^ aijXj X/ . 

For i > n — k, the fth equation of (x', 0/^) o Ax = (x', Ojt) simplifies to = 0, 
which is trivially satisfied. □ 



We will show that we can avoid generating all 2" prefixes. 

Let ('"j^*^') denote the family of vectors (xj, . . . , x„_/f) e {0, 1}"-*^ for which 
Xi + ■ ■ ■ +x„_i- = k. Let denote the family of vectors (xj, . . .,x„_j^) e 

{0, 1}"^*^ for which Xj + ■ ■ ■ + x„_i^_i = k and x„_]^ = 1. The notation is moti- 
vated by the cardinalities. 



[n-k] 
k 



n — k 
k 



and 



■k-1] 
k 



Lemma 3. Let x E {0, 1}". Then there is exactly one k E {0, . . . , [n/2\ } such that 
the {n - k)th prefix x' = (xi, . . . , x„_/,) belongs to {^"~''^) or ('"^^"^')^. 

Proof. Let ivj^ = xi + ■ ■ ■ + x„_j-- Choose the smallest k > with the property 

(i) lUjf = k, or 

(ii) = k + l, x„_/f = 1. 

To see that such k exists consider the pairs of values {wf;, k) for k = 0,1, . . .. The 
first pair is (xj + ■ ■ ■ + x„,0). The left value decreases monotonically by or 

1 down to 0. The right value increases 0, 1,2, The only possibility that (i) 

lacks a solutions is that ivjf and k pass each other, i.e., there exists k such that 

w/^ = k + l,wi^^i = k. 
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Figure 1: The sets ('"^''1) and ('""^"^') ^ for « = 7 and fc e {0, . . . , \n/2\ }. The 
bottom row shows all 128 bit patterns on 7 bits. 



But in that case x„_i- = 1 so that {ii) is satisfied. Note also that if (n) holds for 
k, then (z) does not for that k, and vice versa. Also note that k is at most [n/2j, 
since w^^ is at most n — k. 

We have that either (f) is satisfied and x' belongs to ('''^'^'), or (//) is satisfied 
and x' belongs to ('""^^^1)^. 

We claim that there cannot be another solution k' > k to either (/) or (ii). 
Consider first the case that k' satisfies (z) and k satisfies {ii). In particular, k + 1 = 

Wk = + Xn-k = ^/c-i + 1, SO A: = wjt-i- Then, 

k' = iVjft < W](_i = k, 

a contradiction. The remaining three cases follow from the monotonicity of the 
left and right hand sides. □ 

We note in passing that the previous lemma gives a combinatorial proof of 
the identity 




see figure [T] For another way to verify this expression, compute 

which comes from evaluating a known, closed form for the generating function 
X]jt ('^^'O^*^ at z = 2 [7. exercise 1.2.9.15]. Then add these values for N = n and 
]V = n-ltoarriveat^(2"+i + (-l)"+i+2" + (-l)") = ^(2 + 1)2" =2". 

Finally, we need to bound the number of solutions to (|51l. A priori, this is a 
system in k variables each ranging over {0, 1}, so the number of solutions can 
be 2^. However, we will show that in expectation, the number of solutions is 
exactly 1 : 

Lemma 4. Let x' denote a vector of fixed values x' = {x^, . . . , x„_i^) from ('"^'^') or 
^"'^ introduce the variables x„_t-_|_i, . . . ,x„ E {0, 1}. Setx = (xi . . . , x„). 
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Let A be a matrix whose diagonal entries are chosen uniformly at random from {0, 1}. 
Then the expected number of solutions to 

(x', Ofc) o Ax = (x', 0/c) (mod 2) 

zs 1. 

Proof. Set r, = a,,. For < i < n — k, the ith equation has the form 

n 

Xi ^ij^j = (mod 2) . 

Since there are exactly k entries x, = 1, there are exactly k of such equations 
that do not trivialise to = 0. We can rewrite each of these equations to 

n 

aijXj = 1 (mod 2) 
or, isolating the variables on the left hand side, 

n n—k 

^ij^j = l + ri+ ^ aijXj (mod 2) . 

j=n-k+l j=l 

Now it is clear that the system consists of k equations in the k variables 
. . . , x„. The right hand side is uniformly distributed from {0, 1}, so an assign- 
ment to the variables satisfies an equation with probability ^ • By independence 
of the choices of r,, all k equations are satisfied with probability There are 
2^ assignments, so the expected number of solutions is 2^2^^ = 1. □ 

Algorithm C {Contributing subsets for general graphs.) Given a graph G = {y,E) 
with n = \V\, list the vertex subsets X C V for which dx{X) is odd for all 
xeX. 

CI [Add random self loops.] Choose r^, . . . ,r„ E {0, 1} uniformly at random 
and add a self -loop at vertex if r; = 1. Let A denote the adjacency 
matrix of the resulting graph. 

C2 [Initialize.] Set A: = 0, x' = 0„. 

C3 [Solve.] Solve the nonhomogeneous linear equation ^ in the free variables 
^n-k+l' • ■ • / ^« using algorithm G. If there is no solution, proceed to C5. 

C4 [Filter solution.] Generate every vector x = -|- ai^'i + ■ ■ ■ -|- oi^b^ for all 2'^ 
choices of ai, ■ ■ ■ , G {0, 1}. If x solves the equation x o Ax = 1„, then 
output X. 

C5 [Next prefix.] Update x' and possibly k so that x' is the next element of 

('"t*^') o'" ^ ~ 0/ • ■ • / ^rid return to C3. (When all these 

elements have been generated, the algorithm terminates.) 



11 



Lemma 5. The running time of algorithm C is within a polynomial factor of 



in expectation. 

Proof. For the running time, we first observe that the members of ('"^'^') and 
^[n-k-i]^ can be listed with polynomial (in fact, constant) delay (7| sec. 7.2.1.3]. 

The number of iterations is given by the sizes of (^"^'^') and ('"~^~^^)-^, 
summed over all A; = 0, . . . , [n /2J . These are well-studied quantities, and it 
is known |Zl ex. 1.2.8.16] that 

where F„ is the nth Fibonacci number, bounded by 

The running time is dominated by the time spent in steps C3 and C4. Step 
C3 takes polynomial time every time it is executed thanks to the polynomial 
running time of algorithm G. Step C4 looks more difficult to boimd, because 
it may exhaust the entire solution space to an equation system in k unknowns. 
Naively, this would lead to a total running time of 



n — k 
k 



2k > 



However, by lemma S) the expected size of this solution space is 1. Thus, we 
expect to generate and verify only one vector x, so the time spent in C4 is 
polynomial in expectation in each iteration. By linearity of expectation, the 
total time spent in C4 is within a polynomial factor of l|9]|. □ 



4.1 Derandomization 

We use the method of conditional expectations to derandomise algorithm C. 
Our randomness is over the choices of rj, . . . ,r„, the self-loop indicators. The 
expected number of candidate solutions S given a specific choice of values for 
the first k choices satisfies 

E[S| (ri r,)] = iE[S| (ri,...,r,,0)] + iE[S|(ri,...,r,,l)] 

for all > 0. In particular, there is a specific assignment to r^, . . . , r„ such that 
E[ S I {ri, . . . , r„) ] < E[S]. Thus, by mimicking the behaviour of algorithm C, 
but computing the expected solution sizes of the systems l|8ll instead of solving 
them, we can choose the assignment greedily one variable at the time, until 
we get rid of all the randomness. The final assignment has no more candidate 
solutions than the original expected number. 

The next lemma shows that we can efficiently compute the expectations 
E[S I (ri,...,rfc,0)]andE[S | (ri, . . . ,r^, 1) ]. 
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Lemma 6. Let x' and x be an in lemma^ Let I G {0, . . . , n}. Let A he a matrix 
whose diagonal entries a,, = r,- for I < i < n are chosen uniformly at random from 
{0, 1}. Then the expected number of solutions to 

(x', Ofc) o Ax = (x', 0^) (mod 2) 

is given by a polynomial size linear equation system in the variables ■ ■ ■ ,Xn 

and ri+i,...,r„. 

Proof. The only contributing equations are those for which x, = 1, of the form 

n 

1 • aijXj = 1 (mod 2) . 

We can rewrite these equations as 

ajjXj = ^ + J2 ^iFi (i^od 2) 
j=n-k+l ;=1 

for i < I, and 

n n—k 

ri + Yj ^ij^j ~ 1 + X!( ^v-*-; (n^od 2) 

j=n~k+l j=l 

for f > /. Now it is clear that the system consists of k equations in the k+ {n — I) 
variables x„_j-+i/ ■ • ■ , and r;^x, ...,/'„. 

If this system has S solutions then the expected number of solutions to (|8)l 
isS/2"-'+i. □ 

Algorithm D {Derandomisation of C.) Given a graph G = {V,E) with n = \V\, 
determine values r^, . . . ,r„ G {0, 1} for algorithm C. Let A denote the adja- 
cency matrix of G. 

Dl [hiitialize.] Set / = 1. 

D2 [InitiaUze.] Set A: = 0, x' = 0„, Nq = Ni = 0. 

D3 [Solve for each choice of r;.] For b = 0, tentatively set o;; = b and solve the 
nonhomogeneous linear equation ^ in the free variables . . . ,x,i 

and r;_|_i, . . . , using algorithm G. Let Si, denote the size of the solution 
space. Repeat step D3 for b = L 

D4 [Tally solutions.] Increase Nb by Sb/2"-'+i for b = 0, 1. 

D5 [Next prefix.] Update x' and possibly k so that x' is the next element of 

(["^''1) or ('""1^"^') for fc = 0, . . . , ^n, and return to D3. (When all these 
elements have been generated, proceed to D6.) 

D6 If No > Ni then fix o;; = 1, otherwise fl;; = 0. Increase Z. If / < n return to 
D2. Otherwise output flu, . . . , 
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5 Proof of Theorem |3] 



Let ^ C 2^ denote the set of (directed) cycle covers of the directed graph G. 
We use Iverson's bracket notation: for a proposition P, we write 

[P] = l'' ''''' 

I 0, otherwise . 

G [Y] is the graph induced by Y C V . We sometimes write e G G [Y] to refer 
to a directed edge e in the induced graph G[Y]. We will operate on sums of 
partitions of vertex sets, and manipulate these sums. For this purpose, we 
introduce the notation 

V 

E 

for the sum over all ordered _fC-partitions Yi, . . . , Y^ of V. That is, Yj U ■ ■ ■ U 
Yk = V and Y; n Yy = (1 < / < ; < k). Parts may be empty. 
For a subset T C E of directed edges define 

Mr) = ^ E n[^eG[Yi]u---uG[YK]]. 

Yi,...,YKeeT 

We will show that the residue modulo X of the function h serves as an in- 
dicator variable for Hamiltonicity on This can be viewed as an instantiation 
of the Cut-and-Coimt technique from 13 . 

Lemma 7. For C e 

[Ce^]=/i(C) (modX). (10) 

Proof. First, we show that for C G we have h{C) — 1 (mod K) . Consider a 
Hamiltonian cycle C and a partition Yi, . . . , Y^ of V. First, every partition that 
sets Yj; = y for some k e {1, . . . , X} (and all other parts empty) has e e G[Yj;] 
for all e G C, so its contribution is 1. On the other hand, consider a partition 
where some Yj^ (1 < fc < X) is neither nor V. Since C is Hamiltonian, there is 
a directed edge uv E C with u E Yj. and v ^ Yj.. This directed edge belongs to 
none of the G [ Yj.] , so the product vanishes and the partition does not contribute 
to the sum. Thus, the total contribution of C is K. 

Second, we show for each C i Jif that h{C) is (mod K). In particular, we 
show 

V 

E nt^e G[Yi]U---UG[Yk]] =0 (modX^). 

Yi,...,YKeeC 

Partition the non-Hamiltonian cycle cover C into cycles 

C = Ci U • • • U Cr , 
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such that each C, is a simple cycle. Note that r > 1 because C ^ . 
Construct the corresponding partition of vertices 

V = ViU...UVr, 

such that Vj are the vertices visited by the cycle C, . 
Let Yi, . . . , Yk be a vertex partition such that 

Yl[eEG[Y,]U---UG[YK]] = l. (11) 

eeC 

We first need to observe that the cycle-induced vertex partition ¥■[,...,¥,■ re- 
fines Yi,. . .,yK. Indeed, assume that Vj D Yj is neither empty nor V^. Then 
there is a directed edge e = iiv & Ci with ii e Vj n Yj but v ^ ViCi Yj. In 
particular, the directed edge e belongs to neither G [Yj] nor to any other of the 
G[Yi^], contradicting i fTTIl . Hence the partition Yj, . . . , Yj< consists of the r parts 
of Vi, . . . , y,-- There are K'' ways to pick these parts, each amounting to one 1 in 
the summation. Since r > 1 this shows that the total contribution of C ^ is 
a multiple of and after division of K still vanishes modulo K. □ 

In particular, we can count the number of Hamiltonian cycles modulo K as 

1^1 = J2 HQ (mod K) . (12) 

We next rewrite the right hand side by an application of inclusion-exclusion. 

For a vertex set Z C V let ^(Z) be the family of edge subsets in which 
every vertex has outdegree 1 and all terminals are in Z. In other words, ^(Z) 
is the family of total fimctions f: V ^ Z where vf{v) is a directed edge in the 
graph for all z; e V. 

Lemma 8. 

j:hiC)= E(-l)'^^''' E HC) (modX). 

C£'^ ZCV CGJ^(Z) 

Proof. Consider a cycle cover C G It belongs to ,^(Z) exactly when Z = V, 
so its total contribution to the right hand side is h{C). By the previous lemma, 
the total contribution of such C G is \J^/f\ (mod K). 

Consider now an edge subset C ^ ^ that is not a cycle cover. Let X denote 
the set of vertices appearing as terminals for the directed edges in C. Then C 
belongs to ^{Z) for every Z with X C Z C y. By the principle of inclusion- 
exclusion, the number of Z with X C Z C V is even (using that X C V because 
C is not a cover). In particular there are as many odd Z as even ones. Thus, the 
contributions of aWC cancel. □ 

Recall that du{X) denotes the number of directed edges from u to a vertex 
in X. For pairwise disjoint subsets Y^, . . . , Yj^^ of V whose union includes all of 
Z, let di,(Z; Yj, . . . , Yk) denote the number of directed edges from v io a vertex 
in Z that stay in the same part as v, formally 

d„(Z; Yi, . . . , Yk) = N(Z n Y^) , where e Y^ . 
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Lemma 9. For Z CV, 

E ^(f) = ^ n MZ) E UMZ;Yi,...,Yk) (modX). 
Proof. Expanding h{F) and rearranging, we have 

E = p E E n[eeG[yi]u...uG[YK]] 

= ^ E E meG[Y^]u...uGm 

= ^ E UMZ}Yi Yk). 

Yi,...,rKV€V 

where the last step is based on coimting in two different ways the number of 
ways that every vertex in V can choose another vertex in Z in the same part of 
Yi Yk. 

The next step is to establish 

Yi,...,Ykv€V ^v€V\Z J Vi YKi'eZ 

(13) 

by induction in | y \ Z | . 

For V = Z the claim is vacuous. Otherwise, select a vertex iv e V\Z 
and split the sum on the left hand side into K sums according to which part 
includes the vertex w. 

f; Y\d,{Z;Yx Yj,) = Si + --- + Sk, 

Yi,...,YkdGV 

where 

V\{w} 

Sfc= E Y[dv{Z;Yi YfcUM Yk). 

Yi,...,YKoeV 

In the product, the factor contributed by the term corresponding to u = a; is 

d^{Z; Yi, . . . , Yfc U { w}, . . . , Yk) = d«;(Z n ( Y^ VJ {w})) = du,{Z^Yk) , 
because w does not belong to Z. Thus, 

V\{w} 

Sfc= £ d^(ZnYO n 4(Z;Yi,...,Yk). 
Yi,...,Yk vev\{w} 
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Summing all the S/^ we arrive at 

S, + ---+Sk^ {(j^dUZnYk)) n d-o{Z;Y,,...,YK)\ 

Yi,...,Yi, I V;c=i / -oeV\{w} J 

V\{zu} 

= d,,{Z) Yl dy{Z;Yi,...,YK) 

Yi,...,Yk veV\{w} 
V\{iu} 

= d^{Z) Y. n d„{Z;Yr,...,YK) 

Yi,...,Yt,veV\{w} 

= ( n '^-(Z)) E Y\dv{Z;Y,r-- ,Yj,), 

^w&V\Z ^ Yi,...,Yi^veZ 

where the last identity follows by induction over w, establishing (iSl l. □ 
We are ready to establish theorem |3l namely, 

ijri = i Y (-i)i^\^ifn^-(^\z))nn^y(^/t) (modx). 

^ Z,Yi,...,Y^ \zeZ / ;c=ly6Yt 

Proof of theorem^ Write Z = V\Z. Combining l(T2ll and lemmas [8] and |9] we 
get: 

1-^1 = ^ E E (-1)1^1 n'^-(2)n'^-(z;^i'---'^j^) (modx) 

= ^ E {-l)\^\Yld-o{Z)Yld,{Z;Yi,...,YK) (modX), 

where we just changed the original summand Z to its complement Z. 

Finally, every y G Z belongs to exactly one Y]^, and for this value of k we 
have dy(Z;Yi, . . . ,Yj() = dy{Y^. Thus, 

Yldy{Z;Yr,...,YK) = Y[Y[dy{Yk)- 

yeZ k=ly(^Y^. 

□ 
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